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Abstract 

■ The goal of research is to devise a modification of the perturbative QCD that should be 

regular in the low-energy region and could serve as a practical means for the analysis of data 
below 1 GeV up to the IR limit. 

The proposed "massive analytic pQCD" is constructed on the two grounds. The first is the 
pQCD with only one parameter added, the effective "glueball mass" m g i ~ 1 GeV , serving as an 
\ IR regulator. The second stems out of the ghost-free Analytic Perturbation Theory comprising 

non-power perturbative expansion that makes it compatible with linear integral transformations. 

The differential recurrent relations connecting expansion functions provide a new possibility 
for defining an individual scale for a given physical process. 
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o ; 1 Motivation and Outline 

As it is well known, the so-called perturbative QCD (pQCD) or the renorm-group(RG)- 
improved QCD perturbation expansion taken in the UV limit is a firmly established part 
1 of the current particle interaction theory. This piece is not only respectable but worthy of 

admiration as, starting with gauge-non-invariant quantization, it correlates several dozen 
■ °f experiments at quite different scales from a few up to hundreds of GeV . 

At the same time, the pQCD meets serious troubles in the low energy (large distance) 
domain below a few GeV at the scales marked by the QCD parameter A ~ 400 MeV . This 
Achilles' heel is related to its UV origin. 

To avoid the unwanted singularity in the LE region, several modifications[lJ of the 
pQCD have been proposed. Recently, one of them, the Analytic Perturbation Theory[2J 
(APT), was good enough [5] in describing the polarized Y p l ~ n {Q 2 ) = Ti(Q 2 ) form-factor of 
the Bjorken Sum Rules (BjSR) amplitude down to a few hundred MeV . 

The Bjorken moment is presented there as a sum of PT and higher twist (HT) contri- 
butions 
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+ T HT ; r^T-J^^p^, (1) 



with Agj 7 ^ including the forth-order term ~ (a s (Q 2 )) A . However, an attempt to fit JLab 
data by expression (1) with appropriate HT coefficients failed as the perturbative part 
exploded in the region 0.5 - 1 GeV and the extracted (via comparison with fitted JLab 
data) /i 2i values turned out to be unstable w.r.t. higher loop terms in the first PT sum. 
This prevents the description of data below 1 GeV . 

Along with eq.flT]), in [5] the PT sum was changed for the APT one 

A^ = J2c k «Q 2 )) k =► Kr(Q 2 ) = J2 C *MQ 2 ) (2) 

fc<4 fc<4 
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with Ak{Q 2 ) , the APT ghost-free expansion functions. The positive result consists in good 
fitting of the precise JLab data down to a few hundred MeV with stable HT parameters. 

This achievement gives hope for the possibility of a global fitting down to the IR limit. 
Unhappily, none of the above mentioned ghost-free modifications is suitable for this goal. 
The common drawback is the use of UV logs in the IR region. 

To approach the global fitting of data (like ones for the BjSR form factor), one needs 
to have a theoretical framework with two essential features: 

• Correspondence with common pQCD in the UV (that is above a few GeV); 

• Correlation with lattice simulation results for the effective coupling a s (Q) ghost-free 
behavior in the low-energy domain. 

As a primary launch pad for this construction, the above mentioned APT seems good. 
It satisfies the first of the conditions and, qualitatively, the second one. To exempt the 
APT-like scheme from its last drawback - the singularity (with an infinite derivative) in 
the IR limit, one has to disentangle it from the UV logs. To this goal, a mass- dependent 
RG-invariant modification inspired by our paper [6] will be used. 

In Sect. 3, on the basis of the massive renorm-group (see, Sect. 2), the non-singular 
version of pQCD with one additional (besides A) parameter, an effective "gluonic" mass, 
m g i , a massive pQCD - MPT for short, is formulated. 

2 Massive Perturbation Theory 

In [6], a particular way of constructing the QCD invariant coupling a s (Q 2 ) free of un- 
physical singularities was proposed. In contrary to the APT, it does not involve explicit 
nonperturbative contributions. Instead, the Q 2 algebraic (non-log) dependence appeared 
there due to threshold effects, and an essential technical ingredient was the assumption of 
the finite gluon mass formal presence. 

The model expression for a s (Q 2 ) was obtained there by the RG summation of the 
mass-dependent diagram contribution - see below Eqs. (jlj) and fl6]). It depends upon the 
gluon m g i and light quark mi q (Iq = u, d, s) masses; in the IR region Q 2 > has no 
singularities with a finite limiting value a s (0) and, as Q 2 /m 2 — > oo , smoothly transits 
into the usual asymptotic freedom formula. 

2.1 Mass-dependent 1-loop diagram 

At the one-loop level, the starting element is the massive (mass-dependent) 1-loop contri- 
bution. For example, to the virtual dissociation of a vector particle (photon, gluon) into a 
massive fermion antifermion pair (e + + e~ ; q + q ) in the s-wave state, there corresponds 
a function I s (Q 2 /m 2 ) representable via spectral integrals 

oo 

^) = */lSf; W = ^F; /„(0) = 0.492, />) = 2/3. (3) 
1 

1 For an explicit expression see Sect. 24.1 in the text-book [7] and Sect. 35.1 in the monograph [8]. 
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which in the space-like region z > is a positive, monotonically growing function with 
the log asymptotic behavior I^ v (z = Q 2 /m 2 ) ~ In z — C s + 0(1/ z) ; C s = 2(l — In 2) and 
the regular IR limit with a finite derivative. 



2.2 Massive Renorm-group summation 

For the QCD coupling modification at small space-like Q 2 < A 2 we involve the mass- 
dependent Bogoliubov renorm-group (mRG) formulated in the pioneer RG papers [9] in 
the mid-fifties. As it is known since that time [10], the mRG, like the common massless 
RG, sums iterations of a one-loop contributiorH a s (^)pt = a s — a 2 Ai(z,y) + . . . , into the 
geometric progression 

^(Q 2 )S = n , ""I, r ; z=%, y=4- (4) 



For the 2-loop case with 

a.(z,y) l $ =a s -a 2 s A l {z,y) + al{A\-A 2 {z 1 y)) + ... , (5) 
an analogous (approximate) RG-invariant "massive" sum 

®-s{z,y)f} = as A -, \ • (6) 

hg i + as A 1 (z,y) + a s ^ ) (l + a s A 1 (...)) 

was also devised later (eq.(8) of paper [H]). There, 

A^ 2 \z,y) = l\Q 2 /m 2 )-l\ l ?/m 2 )\ uv ^ In (Q 2 //i 2 ) , (7) 

with transition to the QCD scale , 

« s (Q 2 )ft ] ^ 



fa HQ 2 /A 2 ) 
performed via the relation 

l//3 « s + ln(g 2 //i 2 ) = ln(g 2 /A 2 ). (9) 



3 The MPT construction 

3.1 One-parameter massive Model 

A simple idea is to change the usual UV logarithm \\ix (x = Q 2 /A 2 ) that is singular also 
in the IR for the "long logarithm" L^(x) = ln(£ + x) which reproduces qualitatively the 
smooth LE behavior of function ([3]) being regular at Q 2 = . It is noteworthy that the new 
parameter is expressed via the coupling constant (at Q 2 = 0) by the relation £ = e 1 ^ as 
non-analytic at a s = . On the other hand, it corresponds to the "effective gluonic mass" 
m gi — a/£ A , a rather old notion (see, a fresh survey by Simonov Ref. [12] and references 
therein) used as an effective IR regulator. 

2 Here and below, the superscript in square brackets denotes the order of loop approximation. 
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One-loop case. The "1-loop structure" in the denominator of eg. (151) is changed now 
for the "long logarithm" with £ , an adjustable parameter 

lnx -> L((x) = ln(f + x), x = Q 2 /A 2 . (10) 

At moderate LE scales the form 

= ln£ + ln(l + 0x) = ^_ + l n (l + X ) ; = l/£ (n) 

is more adequate. In terms of this LE form, one has 

aN ( r - c )- 1 - — ■ ch-p- 1/)3oas (12) 

AMPTlX ' 4j "/3 o ln(£ + a:)-l + a s /3oln(l+a;0)' 0-6 ' {U) 

with a s = a s (x = 0); a s |^=io±2 = 0.61 =F 0.05 . The derivative at Q = also could be of 
interest A^\ IPT (0) = —f3 a 2 . 

The 2-loop 1-parameter model. Starting with the 2-loop massive RG-summed result, 
eq.Q corresponding to the mass-dependent PT expansion, eq.floT), we generalize eq. (ll2j) 
by using the same "long logarithm" model for the second-loop contribution 

A 2 {x4>) = m(l + </>x). (13) 

That is 



✓4.1 A/PT 1 (^-') — Q (14) 

' 1 + a s p ln(l + <f> x) + a s § ln[l + a s (3 ln(l + <j> x)] . 

Now, the condition a s {M 2 ) = 0.34 can be used for a rough evaluation of an "effective 
MPT-QCD scale" value. At the NLO case, the A [2] (f) depend ence on £ turns out to be 
rather weak El : 

A [2] (10±2) ~ 315 T- 10MeV. (15) 
Note that the related m g i = \/£ A! 2 ' values could be close to the nucleon mass at £ ~ 10 . 

3.2 Higher MPT expansion functions 

In the construction under devising, we intend to preserve an essential APT feature, namely, 
the non-polynomiality of the modified "perturbative" MPT-expansion, expansion over a set 
of functions {Ak(Q 2 )} connected by differential recursion relations of the type Ak+i(x) :~ 
Ak(x) . To the arguments ascending to the 80s [14] and related to the 7r 2 -terms sum- 
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mation procedure in the s-channel (see, also Refs.[2]) one can add a fresher reasoning 

The "MPT coupling squared". Let us try to use the same recurrent relation (with the 
dotted notation for logarithmic derivative F = x F ) 

- f ^Ak,MPT(x) ■= ~\ Ak,MPT = A) A k +l,MPT(x) + /?1 Ak+2,MPr(x) + ... , (16) 

as in the APT (and logarithmic PT). 



3 with values less than the pQCD one A™'= 3 = 420 ± 10 MeV . 
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This recurrence property ensures compatibility [TB] with linear transformations and, 
particularly, with the transition to the distance picture (Fourier-conjugated with the mo- 
mentum-transfer one) and to the annihilation s-channel. 

In a particular case k = 1, with (fl2|) . and neglecting the second r.h.s. term, i.e. using 
the one-loop relation 

Ak+l,MPT(x) — — Ak,MPT (17) 

as a definition for higher functions, one gets the second expansion function^ 

A 2 , M pt(x) = (A^mpt) 2 R(x) ; R(x) = = J^— , (18) 

o In x 1 + <fi x 

that turns to zero in the IR limit. Besides, A 2 mpt(O) = a s ^ • 

The third MPT expansion function obtained by eqs.( TT7|) and ( f!8|) can be represented in 
the form 

- 2/3 A^ MPT (x) = A 2 ,mpt(x) = 2Ai,mptA 1 , mpt R(x) + A\ MPT R{x) , (19) 
that is sufficient for defining its IR properties A^ MPT (0) = 0; A^\, IPT (0) = —a 2 (p. 

4 General features of the MPT scheme 

Here, we shortly discuss some important features of the proposed MPT construction. 
4.1 Annihilation channel 

For the transition to the s-channel, f{Q 2 ) — > F(s) one uses the spectral representation and 
the Adler-type relation 



f(Q 
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if^, im=Q*r F{s)ia 
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[s + Q 2 ) 2 



which, like in the APT, results in 



1 f°° rlrr 

F(s) = - —p(a); p(a) = Qf{-a-ie). (20) 

Js cr 



In the simplest one-loop case (fT2|) . 

a a [1] (r-n- 1 - 1 r P^ da n ( n \- m , 

1 r°° dn 1 1 ( S TJl 2 \ 

hx&tprtet) = -J —PiW = - aictan TJT) ; Ms) = ln yl^J ■ (22) 

This expression remains regular around s ~ m 2 gl + A 2 and is constant below s = m 2 gl . 
4 For discussion of a more accurate definition of the two-loop second function, sec Appendix. 
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5 Comparison with APT and Discussion 



To compare the new construction with the APT one, in Fig.l we give the curves at a few 
values of £ = 10 ± 2 in the region below 2 GeV for the first MPT function Af\. IPT (114j) 
together with the corresponding APT (dashed) curve. 
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Figure 1: The first MPT function (for a few £ values) in comparison with the APT one. 

It can be seen from the NLO curves that values £ = 8 — 10 seem to be preferable. 
Indeed, for these values, the first Af\, IPT (x) function is reasonably close to the first APT 
one down to 1 GeV . At the same time, in the region around 500-700 MeV it deviates 
from APT but is more similar to the results of lattice simulations, especially to the Orsay 
group [T7l [18] ones. 
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Figure 2: The second and third MPT functions vs. the APT ones. 

Figure 2 exposes the second and third MPT functions as roughly estimated by the one- 
loop Eqs. frT8|) . f|T9|) . It is seen that instead of singular IR slopes of all APT thick dashed 
(red) curves all the MPT functions are the IR regular ones. At the same time, just as in 
APT, the second and third MPT functions A2&MPT are noticeably smaller than square 
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("4i,mpt) 2 and cube (A^mpt) 3 of the first one. E.g., at Q = 0.5 GeV and Q = 1.0 GeV 
the reduction factors for the A^-mpt are about 0.2 and between 0.1 - 0.2 correspondingly. 
Besides, the A^mpt is negative below 1 GeV. 

Now, the "MPT-perturbative" expansion similar to eq.flSJ), due to recurrent relation 
( flTl) can be represented in a form of Taylor series expanded over the parameter A In x = 
Ciir/fto = 1.60 with the final effect 

AW-MwptM; x*=xe- Aln * = Q7(A*) 2 , A* = 2.25 A. (23) 

One should keep in mind that both the logarithm shift A In x and the "one-term approxi- 
mation" ( I2"3"j) error 5A MPT ~ A^,mpt are scheme- dependent quantities. In the MS scheme 
under consideration, SAmpt ~ due to the smallness of the reduction-factor - is negligible. 

However, one can get another angle on Eq. (1231) and return to the old idea of the effective 
coupling constant [IH] which is not so far from "RESIPE" [20] and from the "commensurate 
scale relations" [21] concepts. Then, the new scale can be treated as a specific one for the 
given process; A* = Asj ■ 

Table. "Glueball mass", as estimated for few values of £. 
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729 
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249 
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315 


995 


709 


12 


253 


876 


305 


1160 
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The "glueball mass" values given in the Table for the LO and NLO cases also look 
attractive. They can be confronted with the glueball mass M®® ~ 1 — 2 GeV of paper [12] 
and with gluon mass M ~ 500 MeV from the lattice estimate [17] as well as from solution 
of the Schwinger- Dyson equations (see Ref. [22] and references therein). 

Besides, as it can be shown [23], the MPT perturbative sum Ampt does allow to fit the 
JLab data (together with the duly modified HT contribution) down to the very IR limit. 
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Appendix - 2-loop MPT higher functions 

For a more accurate defining the 2-loop higher expansion functions, one should use 
recurrent relation (|T6|) at k = 1 and truncated eq. (fT7|) for the k = 2 case. 
With temporary technical notations 

4>(t = Inx) = -^A^mpt(x) ; A [ ^ MPT (x) = y(t) , 
one gets two relations of eq. flTB"]) . Neglecting by A^mpt we come to boundary value problem 

y(t)-9y(t) = <P(t), y(oo) = 0; 9 = ^ (24) 
and auxiliary relations 

A^mpt^) = ~Yf > Aa/pt(x) = , (25) 

Solution of ( jUJ) being expanded in powers of 6> 

/"OO 

y(t)= e- s (f)(t + s6)ds = (j>{t) +6(j)(t) + O(# 2 0) (26) 
Jo 

can be rewritten in the form 

J%mpt(?) = -i A?!mpt(x) - % Mmpt + 0(A iMPT ) , (27) 

completely correlating with eq. (116j) . 

In the context of eq. ff23l) . due to the 3-flavour numerical values 0o(nf = 3) = 0.7162 ; 
Pi/Po = 0.5659 , 9 = 0.395 , we note that the second term in the r.h.s. of the last expansion 
- correcting the approximation ()18p - reduces further the error 5 A* of expression (T23]) for 
A* 

MPT ■ 
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